NASA TECHNICAL 

EMORANDUM 




8 

I 1 ? 

X 


< 

co 

< 


GPO PRICE $ . 
CFSTI PRICE(S) $ . 

Hard copy (HC) 
Microfiche (MF) 

ff 653 July 65 



(NASA TM X- 52017^ 


N65-35 243 

C F S SION NUMBER! 


(ACCESSION NUMBER) 

jlZ. 


IASA CR OR TMX OR AD NUMBER) 


(CATEGORY) 


GRADUAL TRANSITION OF NUCLEATE BOILING 
FROM DISCRETE-BUBBLE REGIME 
TO MULTIBUBBLE REGIME 


Is^l^ih-Yun Hsu 3| 

flSf), Lewis Research Cent er ^~T) 
Cleveland, Ohio J 





Ph€ |j>rir)f 


v TEUH N LOAl, PKEPKUU piepared -f»r N Heat Transfer 
Confer enee and Prophets Show sponsored by the. 

Am erican Inst^teW of Chem freal E ngjf jie c * o and 4j§e 
Am e r ican Societ y of Me ch anical Engirtoera 
Cleveland, Ohio,rAuS^T 9-12, 1964 ' 

y 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D.C. * 1964 


TECHNICAL MEMORANDUM 


GRADUAL TRANSITION OF NUCLEATE BOILING 

FROM DISCRETE-BUBBLE REGIME 

TO MULTIBUBBLE REGIME 

by Yih-Yun Hsu 
Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PREPRINT prepared for 

Heat Transfer Conference and Products Show 
sponsored by the American Institute of Chemical Engineers 
and the American Society of Mechanical Engineers 
Cleveland, Ohio, August 9-12, 1964 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


GRADUAL TRANSITION OF NUCLEATE BOILING 


FROM DISCRETE-BUBBLE REGIME 
TO MULTIBUBBLE REGIME 
by Yih-Yun Hsu 
Lewis Research Center 

ABSTRACT ^^3 

A photographic study was made of over 5000 bubbles in the nucleate 
boiling of methanol and water on a narrow heating strip at various heat 
fluxes and degrees of subcooling. The result showed that transition from 
the discrete-bubble regime to the multibubble regime was gradual. The 
fraction of heating area covered by multibubbles increases with the in- 
creasing heat flux and is predictable. The area fraction is a Poisson 
function of the product of the mean area of i 
and the instantaneous population density. 
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INTRODUCTION 

Since the nucleate boiling is recognized as a very effective means 
of heat transfer, a tremendous amount of effort has been directed toward 
the understanding of this interesting phenomenon. Among such efforts, 
a good part has been devoted to the study of bubbles. To facilitate 
observation, studies usually were made on discrete bubbles. Thus, the 
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formulation of theories on nucleate boiling were based upon the informa- 
tion derived from discrete bubbles (refs. 1 to 3). These theories were 
generally applied to the entire regime of nucleate boiling. However, 
recently it has been becoming more and more clear that there actually 
exist several subdivisions in nuclear boiling; namely, a discrete bubble 
regime, a merging bubble regime, and perhaps a vapor-patch regime (refs. 4 
to 7). It is evident that theories should be developed to deal with each 
region individually, as well as to predict the transition from one region 
to another. 

It is the purpose of this paper to study the transition of the 
discrete-bubble regime to the merging-bubble regime. In reference 8, an 
abrupt transition point was proposed, while in reality the transition is 
gradual and continuous. This paper will show how the area covered by 
merging bubbles gradually increases with increasing heat flux, and that 
the area fractions for merging bubbles can be related with other param- 
eters such as bubble size and instantaneous bubble population. The hope- 
is that if the area fractions covered by the merging bubbles and the 
discrete-bubbles at a given condition are known, the overall heat- 
transfer coefficient can be synthesized by weighting the contributions 
due to the two bubbling mechanisms according to their respective area 
fractions. 

The experimental phase of this work consisted of a photographic 
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study of nucleate boiling of methanol and water on a — - by --inch heat- 
ing strip under 1 atmosphere pressure. The resulting data were then 
analyzed by assuming a Poisson distribution of bubbles. 
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NOMENCLATURE 

A area 

a empirical area parameter used in equation (l) 

a,b empirical parameters used for bubble growth rate R = at b 
D bubble diameter 

f bubble generating frequency 

g gravitational acceleration 

h total number of all single bubbles studied in one roll of film 

K thermal conductivity 

k number of sample frames 

L length of heating strip 

M average number of sites per cell 

N site population 

n instantaneous bubble population 

P probability according to Poisson function 

q heat flux 

R bubble radius 

R bubble growth rate, dR/dt 

s standard deviation associated with average area fraction of merging 

bubbles 

AT su b subcooling temperature difference between saturation and bulk 
temperatures 

ATy temperature difference between surface and bulk 
t time 

t K bubble growth period 
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t w waiting period 

W half-width of heating strip 

X number of bubble sites per cell 

x number of bubbles in a cell 

g contact angle, radians 

S thermal layer thickness 

A latent heat of vaporization 

p average number of bubbles per cell 

p density 

a surface tension 

cp area fraction 

Subscript: 
av average 

b bubble base 

B bubble 

calc calculated 

d departure 

exp experimental 

F Fritz equation of bubble departure, equation (9) 

Z liquid 

m merging bubbles or multibubbles 

s single bubble 

sub subcooling 

t total 

v vapor 
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Literature Survey 

As mentioned in the INTRODUCTION* bubble interference has been re- 
ported previously. The earliest mention of it was probably that found 
in reference 9. More recent experimental findings have been reported 
since then. In general* the bubble interference can be classified into 
two types* those of vertical interference and those of lateral interfer- 
ence. The vertical interference occurs between consecutive bubbles 
emitted from the same nucleation site in rapid succession. This type of 
bubble interference was called chain-bubbles interference in reference 9. 

It was also reported in references 6* 8* 10* and 11. This type of bubble 
coalescence was the model utilized in reference 8 to derive the criterion 
for the transition from the discrete-bubble regime to the merging-bubble 
regime. Deissler used a similar model for his analysis of burnout heat 
flux (ref. 12). The lateral type of bubble coalescence (or mushroom 
bubbles according to ref. 4) was the interference between the neighbor- 
ing bubbles due to close proximity* a phenomenon not entirely unfamiliar 
to the urban or suburban dwellers. As observed in references 3* 4* and 
13* a growing bubble, while still attached to the heating surface, would 
merge with a neighboring bubble. This merging could be caused either by 
contact of two growing bubbles or by the "up draught" of a departing bubble. 
The area of influence of each bubble has been found to be roughly 2 bubble 
radii away from the nucleation center (refs. 3 to 5). In either case, the 
lateral-merging bubbles could be pictured as mushrooms with two or more 
stems. These stems would be the places where vaporization occurs. This 
type of coalescence has been included in those boiling models postulated 



in references 3 to 5. As will be shown later in the section RESULTS AND 


DISCUSSION, the mushroom bubbles are far more frequently observed than 
the chain bubbles. Therefore, the mushroom bubbles will be the ones dis- 
cussed in this paper. 

Since the lateral coalescence is due to the interference of neighbors, 
the distribution of bubbles should be known. In reference 14, it was 
found from the distribution of the sites on a boiling surface that the 
site population was distributed according to Poisson's equation 

%(x) = e '“( M ) ]C (1) 

A • 

where M = Na/A and X = Na/A, but no attempt was made to predict the 
population distribution a priori; instead, the cell area, a, was used 
strictly as an empirical parameter to fit the data with Poisson curves. 

Apparatus and Procedure 

The test was carried on inside a 6-inch-diameter by 4-inch-high 
cylindrical tank made of stainless steel and provided with viewing win- 
dows. The tank had provisions for a fill, a drain, a pressure gage, 
electrical connections, thermocouple leads, and auxiliary heaters. The 
heater was a thin, electrically conductive, transparent coating l/l6-inch 
wide, 1-inch long deposited on a 1- by 1- by l/8-inch heat-resisting glass 
plate. The plate was mounted horizontally on a small Bakelite bench with 
a mirror situated beneath the plate and inclined 45°. Thus, the camera 
aiming from a front window saw simultaneously a front view and a bottom 
view (through the mirror) of the image of any bubble generated on the 
heating surface. The plate was clamped down by two copper clamps which 
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also served as electrical leads. The actual heating area was l/l6 by 
3/4 inch, since the two end areas were covered by the copper clamps 
(fig. l). The l/l6-inch-wide strip was used instead of a wider heating 
area to assure that no more than two rows of bubbles were generated. This 
arrangement was necessary to avoid confusion in the front view because of 
the presence of overlapping rows of bubbles. Originally, an alternating- 
current source was used, but due to the low heat capacity of the heating 
film, there were 120-cycle-per-second temperature fluctuations on the 
heating surface. At lower heat flux, when bubble frequencies were low, 
this 120-cycle-per-second fluctuation apparently did not have a serious 
effect. However, as the heat flux was increased and bubble frequency fell 
in the vicinity of 120 cycles per second, the alternating- current fluctua- 
tion began to dictate the bubble frequency, and the bubbles began to grow 
in unison. Thus, only those runs where there was no apparent synchroniza- 
tion between bubbles and current waves were retained. Later, a direct- 
current source was used. However, because of the high voltage applied to 
the electrodes (60 to 120 v) , electrolysis would take place if water was 
used. Thus, only methanol was used for direct- current runs. At the be- 
ginning of a series of runs, the tank was loaded with a fresh batch of 
pure methanol or distilled water. The liquid was preheated to a desired 
temperature by the auxiliary heater. The bulk temperature was constantly 
monitored through thermocouple readings, and the temperature level was 
controlled by turning the auxiliary heaters on and off. The auxiliary 
heaters were always off while actual test runs were being carried out. 

The test heater would be turned on, set at a desired heat flux by varying 
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the applied voltage, and then high-speed motion pictures (up to 5000 
frames/sec) vere taken. The simultaneous viewing of bubble activities 
from front and bottom (mirror image) was deemed necessary to get the 
true picture of bubble interaction. However, because of the difference 
in optical distance of the two views, it was extremely difficult to keep 
both views in sharp focus. Thus, the optical qualities of pictures were 
somewhat sacrificed. Because of such difficulties, studies were limited 
to those runs with a moderate amount of merging to maintain the accuracy 
of reading. The high-speed motion pictures were analyzed on a motion- 
picture analyzer. A total of 14 rolls were examined. For each roll, 

50 to 100 frames were studied. The sample frames were selected by arbi- 
trarily stopping the film 50 to 100 times at irregular intervals. 

On each frame, information about each of the bubbles present on the 
entire heat surface 3/4 by l/l6 inch were measured and recorded. The raw 
data include: 

(1) The location and size of each bubble (the size or diameter of a 
bubble is defined as the width of a bubble at its widest part. In a few 
instances, a mushroom bubble could have a very wide hovering bubble over- 
casting a large area. In such cases, the bubble size is defined as the 
width of the bubble stem below the height of an average single bubble. 

See fig. 2. ) 

(2) The classification of bubbles, namely, whether the bubble was a 
single bubble, or was merging with other bubbles. The criterion for the 
merging bubbles was the physical contact of two or more bubbles, while at 
least one bubble was still attached to the heating surface. 
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(5) The number of active bubbles involved in a merging bubble, and 
(4) The total number of active bubbles n on that frame. 

From the aforementioned raw data, calculations were made, and the follow- 
ing was deduced: 

(1) The average size of single bubbles was made over all the single 
bubbles recorded in those sample frames. This was expressed in terms of 
the area fraction of influence of single bubbles, against the total area 
or cp g . The area of influence of the bubble cp s was computed by the 
equations 

*D 2 

cp s = ^ if D s < W (2a) 

and 

2 W-VdI - W 2 + sin -1 t?- 

cp s = — t if D s > W (2b) 

A t 

The shaded areas in figure 3 are the areas of influence of single bubbles. 
Note that equation (2b) represents the area of a part of a circle with 
two segments cut off. The mean area fraction averaged over all the 
bubbles is 

h 

5 Vi 

^s,av = “ H (^c) 

where h is the total number of all the single bubbles studied in the 
roll. 

(2) The average size of the area covered by a merging bubble in a 
frame was expressed as the area fraction cp m , which was calculated by 
summing up all the area fractions covered by each merging bubble in the 
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same frame. The area fraction of each merging bubble was also computed 
from equations (2a) and (2b), except that D m /2 was used instead of D s . 
The average area fraction of merging bubbles was computed over all the 
frames , or 

k 

X ^m,i 

i=l 

™m,av “ ^ 

where k is number of sample frames and cp m ^ is the sum of all the 
area fractions of merging bubbles in the frame i. Also computed for 
each roll of film was the standard deviation s associated with the 

^m^av 



(3) The average instantaneous bubble population density n av was 
taken for the total number of sample frames, k. 

The total number of active sites N seen in a given roll of film 
were also studied. The movie was projected on a paper and all the sites 
where bubbles had ever been generated were marked down. The range of 
conditions and the data are tabulated in table I. 

RESULTS AND DISCUSSION 

General Description of Photographical Observation 
Before the quantitative study on bubble interferences, the following 
qualitative descriptions should be given: 

(l) The merging of bubbles was predominately due to lateral coales- 
cence. The merging took place when one growing bubble got into the area 
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of influence of a neighboring bubble. A detailed listing of bubble 
classification and raw data is contained in reference 15. 

(2) The number of merging bubbles and the area covered by these 
merging bubbles increased with both the increase of bubble size and the 
instantaneous bubble population. With pressure and degrees of subcooling 
held constant, both the bubble size and the instantaneous bubble popula- 
tion increased with heat flux. Thus increasing the heat flux means in- 
creasing the area of merging bubbles. 

(3) The location of merging bubbles appeared to be random. For a 
given location, at one moment there could be anything ranging from no 
bubbles, one bubble, several discrete bubbles, to merging bubbles. How- 
ever, the probability of having merging bubbles increased with increas- 
ing heat flux. 

Analysis 

Based upon the general qualitative description of bubble interfer- 
ence observed photographically, a model will be postulated to account 
quantitatively for the area fraction covered by the merging bubbles. The 
analysis will be carried out in two steps: The first step will be to 

seek the relation between the area of a merging bubble with the quanti- 
ties such as the mean area of influence of a bubble and the instantaneous 
bubble population, provided the latter two are given. The second step 
will be to attempt to estimate the area of influence and the instanta- 
neous bubble population from the more basic information such as heat 
flux, subcooling, and total bubble population (or site population). The 
purpose is to estimate the area of merging bubbles from the aforementioned 


basic information. 
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(1) Basic Model - The basic model for bubble interference is 
described in the following manner: 

(a) Each bubble has an area of influence, which is the area within 
1-bubble-diameter distance from the nucleation center. This assumption 
is based upon the observations made in reference 3. 

(b) Since each bubble grows, the area of influence is based upon the 
time-mean bubble size. 

(c) Two bubbles will merge if one is located within the area of in- 
fluence of the other. 

(d) The bubbles are assumed to be distributed according to a Poisson 
distribution. This distribution is assumed to apply not only to site 
population as found in reference 14, but also to the instantaneous bubble 
population. Note that at any moment, only part of the sites are actively 
occupied by bubbles, while the rest are in the waiting period. 

(2) Poisson Distribution - The following equation is used to express 
the previously postulated model in mathematical form: 

V*) = < 3 > 

K xl 


where P (x) i- s the percentage of cells each of which has x bubbles in 

r 1 

it, while the average number of bubbles per cell is 



n av 


A s = n a v £ Ps 


(4) 


* 


where A s is the mean area of influence of a single bubble or the area of 
a cell, n av is the average instantaneous bubble population on a total 
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heating area A-^, and cp s is the area fraction A s /At> Note that equa- 
tions (l) and (3) are similar in form except that the site population N 
was used in equation (l) , while the instantaneous bubble population n 
was used in equation ( 3 ) and that in equation ( 3 ) the cell area was de- 
fined as the area of influence of a bubble. 

Since it is assumed that bubbles will merge when two or more bubbles 
are present in one cell, the percentage of cells that contain merging 
bubbles is 


P (x 
P 


> 2 ) 




V x) 


(5a) 


or 

P 4 (x > 2) = 1 - P M (0) - P^l) = 1 - e-’WPsU + n av <p s ) (5b) 

Since P )i (x > 2) is by definition the percentage of cells covered by two 
or more bubbles or merging bubbles, it is the area fraction covered by 
merging bubbles 

<p m = P |J (x > 2) = 1 - e-’WPsd + n av cp s ) (6) 


This equation will give the area fraction of merging bubbles if the mean 
instantaneous bubble population n av and mean area of influence of a 
single bubble are known. These two terms can either be obtained experi- 
mentally or analytically. The next two sections constitute the second 
step of analysis, namely, determination of n av and A s analytically. 

(3) Mean Area of Influence of a Single Bubble - According to the 
assumed basic model, the mean area of influence is the area within 
1 bubble diameter from the nucleation center, and the bubble diameter 
is the time average of a growing bubble 
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av 



D(t)dt 


( 7 ) 


or 

i p tg 

R av = T~ / R(t)dt (7a) 

t g Jo 

The term bubble radius R(t) can be obtained through bubble growth 
information. Although many theoretical equations are available, it is 
more convenient to use the empirical expression R = at b , where b = 0. 4 
(ref. 16) will be used. Since the process of computing the time-average 
radius R av involves integration of R(t), small deviations in R(t) 
usually will be evened out. Thus 



t^dt - 


atb 

g_ 

b + 1 


R d 

1 + b 


% 

1.4 


if b = 0. 4 (7b) 


As to the departure radius R^, Staniszewski ' s empirical expression will 
be used 

R d = R f (1 + 10. 44 R d ) (8) 

where Rp is the departure radius according to Fritz* equation 

s *--°- iZ15 *4^r^ t (9) 

is in feet per second, P is the contact angle in radians, and R is 
in feet. 

Unfortunately the growth rate at departure R^ involved in the above 
equation can no longer be calculated from the expression R = at , partly 
because the exponent b actually varies with time and partly because the 
coefficient a should be a function of an experimental condition such as 
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heat flux, subcooling, pressure, cavity size, etc. Thus an expression 
for the growth rate R as a function of the test condition should be 
used. 

Although many bubble growth equations are available, only a few con- 
sidered the effect of the bulk turbulence by including terms that describe 
the thermal layer or the heat dissipation to the bulk. Among such equa- 
tions are those proposed in references 3, 17, and 18. The equation in 
reference 17 would be quite convenient to use if both AT W and q w were 
known. Unfortunately, the growth expressions in references 3 and 18 are 
rather clumsy to use. However, if only the bubble growth rate at depar- 
ture is of interest, the situation is somewhat simpler because of the 
fact that in the later stage of bubble growth the sensible heat stored 
in the superheated layer enveloping the bubble should have been long since 
exhausted. Therefore the bubble should be receiving heat only from the 
bubble base. (This heat may be in the form of evaporation of microlayer as 
shown in ref. 19. ) While losing heat to the surrounding bulk, the bubble 
growth rate at departure can be easily derived by following the procedure 
in reference 3 




( 10 ) 


where — — 0. 25 if the bubble at departure can be assumed as a truncated 
a B 

sphere with contact angle between 45° to 60° (c. f. , ref. 3) and the S is 


the thermal layer thickness. The information about thermal layer thickness 
of a boiling fluid is very meager, but there are a few measurements (ref. 6 
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and 20). Therefore, if the thermal layer thickness 8 is known, by 
using equations (8) to (10), the bubble departure size can be determined. 

(4) Instantaneous Bubble Population - The instantaneous bubble popu- 
lation should be differentiated from the commonly used term "bubble popu- 
lation". The latter actually is a misnomer. When the number of bubble 
columns are counted during the running of a motion picture or the number 
of aureole left on a plate after boiling, only the population of the bubble 
nucleation sites is shown, but not the bubble population at any moment. 

The relation between site population N and the instantaneous bubble popu- 
lation n can be likened to that between the number of houses in a block 
and the number of families actually at home at a given moment. It is easy 
to see that these two populations can be tied by the equation 


n, 


_ g,av 

av “ t + t 

°g,av T °w,av 


N 


- ^av^av^ 


( 11 ) 


The variation of N as function of q has been reported in many places, 
and, unfortunately, the result varies widely. The difficulty stems from 
the diversity of surface condition and hysteresis (refs. 2, 6, and 21). 
Unless some characteristic parameter other than rms roughness of a sur- 
face can be found to account for size distribution of cavity, it is futile 
to try to correlate N against q. However, the site population N is 
still a quantity much easier to determine experimentally than the instan- 
taneous bubble population n. Thus, it is still worthwhile to obtain n 
through N. 

The mean frequency f av can easily be determined through the expres- 


sion in reference 17 
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fD = 0.59 


(p i ~ P v )go 
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so 


f 

oir — 


0.59 


l(p i ~ P v )so 1 / 4 


av ~ n 


d^av 


(13) 


The time of growth period t can readily be calculated through 


R d = 


J g 


R dt 


(14) 


provided R(t) and R d are known. As mentioned before, if an empirical 
expression is used 

(15) 


R d = atg 


R d = abtg -1 


Hd = 

R d b 


(16) 


or 



(17) 


Thus from equations (8), (10), and (17 ) the growth period t can be cal- 
culated. Strictly speaking, equation (10) can be used to replace equa- 
tion (16), only when the empirical form (eq. (16)) is identical to the 
analytical form (eq. (10)). However, an underestimated growth rate R 
tends to give an overestimated growth period t g and an underestimated 
departure diameter D d . The result is that the two errors tend to com- 
pensate each other in the product of mean bubble population and mean area 


of influence 
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n av^s = nD d (18 ) 

(5) Calculation of Area Fraction Covered by Influence of a Single 
Bubble (9 s ,av) ca i c from ^ D av) calc " The mean influence area fraction 
covered by a single bubble can be computed from the time-averaged diameter 
of a bubble D av by using equations (2a) and (2b), except that D av will 
be used in the place of D. The cp s thus computed will be (9 s ,av) C alc* 
Strictly speaking, (9 Sjav ) calc should be an average of (9 s ) calc which > 
in turn, should be computed from D(t) as shown in equation (2c). Since 
only an estimation was intended, ( c P s ,av^ C alc Gan be directly com P ute<i 
from (Dav) calc - 

Comparison of Experimental Data with Analysis 
The quantitative result will be compared with the model derived in 
the section Analysis in two steps also. The first step will be to check 
whether the area fraction of merging bubbles cp m based upon the Poisson 
distribution (eq. (6)) can be used to relate cp m with the experimental 
values of the mean area of influence of a single bubble 9 s?ex p an< ^ 
mean instantaneous bubble population n a v,exp* '^ ie secon( l step will be 
to test whether equations (8) and (ll) can be used to predict 9s, exp 
and n av?ex p respectively and whether the calculated product (9s n av) ca p c 
can be used to predict the merging bubble area fraction 9m, exp’ 

(l) The Relation Between Area Fraction of Merging Bubbles 
Against Product of Measured Values of Mean Area of Influence of Single 
Bubbles and Mean Instantaneous Bubble Population ( n av9 s ) eX p ( i- e * > the 
average number of bubbles per cell) - To test equation (6), (9m,av^ e xp 
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• was plotted against (n av cp s ^ av ) eX p ^^-Sure 4(a). The solid curve repre- 
sents equation (6). Each circle represents the mean values obtained from 
one roll of film. Figure 4(b) shows a similar plot on rectangular coordi- 
nates to show the standard deviations associated with each area fraction 
of merging bubbles for the samples studied. Judging from the figures, the 
model is fairly close. Thus, if the mean bubble population n av and the 
mean area of influence of single bubbles cp g are given, the area covered 
by merging bubbles can be calculated. The values of n av and 9 s ^ av can 
either be obtained experimentally in the same way that figure 4 was con- 
structed, or they can be estimated from test conditions through bubble 
departure size, bubble growth rate, and frequency by using the available 
equations. 

(2) Comparison of Calculated and Measured Bubble Departure Diameters 
D^cai Against D e ,exp - To estimate cp s ^ av , it is necessary to know the 
departure diameter D^. Equations (8) and (9) were used to compute D^ 
with departure growth rate computed from equation (10). The thermal 

layer thickness 6 used in equation (10) was determined from the experi- 
mental measurement in reference 10 by matching the heat- transfer coeffi- 
cient for the case of water (varying roughly around 2X10" ^ ft). For the 
case of methanol, the thermal layer thickness was assumed to be 2X10 -3 
feet. The calculated values of D^ are then compared with the experi- 
mental ones derived through equation (7). The comparison is shown in 
figure 5. It can be seen that most points are within a ±20 percent error 
limit, which is about the same as the 25 percent error limit of equa- 
tion (8). 
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(3) Comparison of Calculated and Measured Mean Instantaneous Bubble 
Population n av ^ calc Against n aV;exp - To test equation (ll), the mean 
instantaneous bubble population n av was calculated from the experimen- 
tally determined site population N together with the calculated fre- 
quency f av and growth period t g ^ av . The bubble frequency f av was 
calculated from equation (13) by using the calculated departure diameter 
®d,cal* The growth period t g ^ av was calculated from equation (17) by 
using Dd. ; calc and Pd,calc (from eq. (10)). The comparison with 
n av,exp shown in figure 6. The +60 percent error limits are also 
shown in the figure, which are the error limits associated with equa- 
tion (12). 

(4) Comparison of Measured Merging Bubble Area Fraction (cp m ) exp 
with That Obtained From the Calculated Average Instantaneous Bubbles 
Per Cell (n av <p s ) calc - By using equations (2a) and (2b), cp s ^ av was com “ 

puted from D a v,calc- In fi S ure 7 > the product (cp s ,av n av) ca i c was 
plotted with the average area fraction of merging bubbles cp m)av - Also 
shown in figure 7 is the theoretical curve from equation (6). Although 
there is scattering, the result is still quite gratifying considering 
all the crude assumptions being made and the large error limits associ- 
ated with empiricism of equations (8) and (12). Thus, it is shown that 
an estimate of the area covered by merging bubbles can be based upon the 
test conditions (heat flux, pressure, subcooling, etc. ) provided that 
the site population is known. 
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CONCLUDING REMARKS 

In the nucleate-boiling regime, if both discrete bubbles and merging 
bubbles are present, the overall heat-transfer coefficient ■will have to 
be determined by considering the contribution due to both bubbling 
mechanisms. Even if the heat-transfer process of each mechanism was 
known, a weighting factor is needed to determine the relative contribu- 
tion of each mechanism. 

One possible weighting factor would be the area fractions of merging 
bubble and discrete bubbles. Based upon the result obtained from the 
boiling of methanol and water on a narrow heating strip, it is found that 
the transition from the discrete-bubble regime to the merging-bubble 
regime is gradual. Furthermore, the results showed that area fraction 
of merging bubble can be predicted satisfactorily from the Poisson dis- 
tribution if the average number of bubbles per cell (the cell is defined 
as the average area of influence of a single bubble) is known. 

It is desirable to be able to determine the average number of bubbles 
per cell a priori. A method of estimating this item based upon crude 
assumptions and empirical equations was proposed in this report. The 
maximum error associated with the estimated values was roughly 100 per- 
cent, which might be due to the large errors introduced in the basic 
empirical equations. If a better method is available for estimating 
the bubble size and bubble population, more accurate estimations of 
the number of bubbles per cell and, thus, of the area fraction of merg- 
ing bubbles might be possible. 
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Figure 1. - Setup of heating element. 
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Figure 2. - Various configurations of bubbl 
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( ) Area covered by bubble with radius R s . 

Y/\ Area covered by influence of bubble ' 

(radius of twice the bubble radius; 
area fraction of bubble influence is <£ s ). 
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(b) Rectangular plot showing standard deviation. 

Figure 4. - Concluded. Area fraction of merging bubble as a func 
tion of measured average number of bubbles per cell. 



Figure 5. - Comparison of calculated and measured bubble 
departure diameters. 
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